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We suggest and demonstrate a scheme to reconstruct the symmetric two-mode squeezed thermal 
states of spectral sideband modes from an optical parametric oscillator. The method is based on 
a single homodyne detector and active stabilization of the cavity. The measurement scheme have 
been successfully tested on different two-mode squeezed thermal states, ranging from uncorrelated 
coherent states to entangled states. 
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Introduction - Homodyne detection (HD) is an effec¬ 
tive tool to characterize the quantum state of light in 
either the time [T]-[8] or the frequency [9ll28] domain. In 
a spectral homodyne detector, the signal under investi¬ 
gation interferes at a balanced beam splitter with a lo¬ 
cal oscillator (LO) with frequency uoq. The two outputs 
undergo a photodetection process and their photocur¬ 
rents are combined leading to a photocurrent continu¬ 
ously varying in time. The information about the spec¬ 
tral field modes at frequencies cjo ± (sidebands) is then 
retrieved by electronically mixing the photocurrent with 
a reference signal with frequency Ft and phase Upon 
varying the phase 0 of the LO, we may access different 
field quadratures, whereas the phase can be adjusted 
to select the symmetric S or antisymmetric A balanced 
combinations of the upper and lower sideband modes. 

Measuring the sole modes S and A through homodyne 
detection is not enough to assess the spectral correlation 
between the modes under investigation [T] and, in turn, 
to fully characterize a generic quantum state. In order to 
retrieve the full information about the sidebands it has 
been suggested that one should spatially separate the two 
modes isni EH or implement more sophisticated setups 
mm involving resonator detection. On the other hand, 
it would be desirable to have schemes, which do not re¬ 
quire structural modifications of the experimental setup. 
In turn, this would make possible to embed more easily 
diagnostic tools in interferometry [33] and continuous- 
variable-based quantum technology. Remarkably, in the 
relevant cases of interest for continuous variable quan¬ 
tum information, such as squeezed state generation by 
spontaneous parametric down-conversion, the correlation 
between the modes vanishes, due to symmetric nature of 
the generated state. 

In this Letter we suggest and demonstrate a mea¬ 
surement scheme where the relevant information for 
the quantum state reconstruction of symmetric spectral 
modes is obtained by using a single homodyne detector 
with active stabilization through the Pound-Drever-Hall 
(PDH) technique |3|. In particular, the reconstruction is 
achieved by exploiting the phase coherence of the setup, 
guaranteed in every step of the experiment, and two aux¬ 
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FIG. I: Schematic diagram of the experimental setup. See 
the main text for details. 


iliary combinations of the sideband modes selected by 
setting the mixer phase at = ±7r/4. 

Homodyne detection and state reconstruction - A 
schematic diagram of our apparatus is sketched in 
Fig-lU The principal radiation source is provided by 
a home made NdiYAG Laser (^300 mW @1064 nm and 
532 nm) internally-frequency-doubled by a periodically 
poled MgOiLiNbOs (PPLN in Fig. [^. To obtain the 
single mode operation, a light diode is placed inside the 
laser cavity. One laser output (@532 nm) pumps the 
MgOiLiNbOs crystal of the optical parametric oscillator 
(OPO) whereas the other output (@1064 nm) is sent to 
a polarizing beam splitter (PBS) to generate the local 
oscillator (LO) and the seed for the OPO. The power 
of the LO (^10 mW) is set by an amplitude modulator 
(AM). Two phase modulators (PMa and PMb in Fig. 
generate both the sidebands used as OPO coherent seeds 
and as active stabilization of the OPO cavity with the 
PDH technique [3]. For the OPO stabilization we use a 
frequency of 110 MHz (HF) while the frequency Fl for 
the generation of the input seed is about 3 MHz. This is 
indeed a major effort, but it will turn out to be funda¬ 
mental for the full reconstruction of the symmetric states 
addressed below. The OPO cavity is linear with a free 
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spectral range (FSR) of 3300 MHz, the output mirror 
has a reflectivity of 92% and the rear mirror of 99%. The 
linewidth is about 55 MHz, thus the OPO stabilization 
frequency HF is well above the OPO linewidth while the 
frequency Q is well inside. In order to actively control 
the length of the OPO cavity its rear mirror is connected 
to a piezo that is controlled by the signal error of the 
PDH apparatus. 

The detector consists of a 50:50 beam splitter, two low 
noise detectors and a differential amplifier based on a 
LMH6624 operational amplifier. The interferometer visi¬ 
bility is about 95%. We remove the low frequency signal 
through a high-pass filter @500 kHz and then the sig¬ 
nal is sent to the demodulation stage. To extract the 
information about the signal at frequency Q we use an 
electronic setup consisting in a phase shifter, a mixer (0 
in Fig. and a low-pass filter @300 kHz. Since, as we 
will see in the following, we need to measure the signal at 
two different orthogonal phases, and T 2 = + 7r/2, 

for the sake of simplicity we implemented a double elec¬ 
tronic setup to observe at the same time the outputs (see 
Fig. 0 . Finally, the LO phase 0 is scanned between 0 and 
27r by a piezo connected with a mirror before the beam 
splitter of the HD. The acquisition time is 20 ms and we 
collect about 100 000 points by a 2 GHz oscilloscope. 

If ao(ci;o) is the photon annihilation operator of the 
signal mode at the input of the HD, it is easy to show 
that the detected photocurrent can be written as (note 
that the “fast term” ujq is canceled by the presence of the 
LO at the same frequency) I{t) oc ao(t)e“'^^ + aQ(t)e*^ 
[ 35 ] . where 0 is the phase difference between signal and 
LO and we introduced the time-dependent field operator 
ao(t), that is slowly varying with respect to the carrier 
at cjo, such that ao(t) = f doo F{uj) ao{uj) = 

e-2u;otQ^o(t), F{uj) being the apparatus spectral response 
function. 

To retrieve the information about the sidebands at 
frequencies cjq i ^5 described by the time-dependent 
field operators d±o(t), we use electronic mixers set 
at the frequency Q with phase shift T with re¬ 
spect to the signal, leading to the current /^(t, T) = 
+ T). Neglecting the terms proportional 
to exp(± 2 ir^t) (low-pass filter), we find the following 
expression for operator describing the (spectral) pho¬ 
tocurrent /o(t, T) oc X 0 (t, T|f 2 ), where ^^(t, T|f]) = 

is the quadrature operator 
associated with the field operator (note the dependence 
on the two sidebands): 




a+a{t)e^^ +a.n{t)e-^^ 

V2 


( 1 ) 


Note that [6(t, ^'1$^)] = X(Ac^)-i (1^ - ^'D- 

The interaction inside the OPO is bilinear and in¬ 
volves the sideband modes a±Q [35]. It is described by 
the effective Hamiltonian Hq oc + h.c., that is 

a two-mode squeezing interaction. Due to the linear¬ 
ity of Hq, if the initial state is a coherent state or the 


vacuum, the generated two-mode state qq is a Gaus¬ 
sian state, namely, a state described by Gaussian Wigner 
functions and, thus, fully characterized by its covariance 
matrix (GM) ctq and first moment vector R [2l|37j- It 
is worth noting that due to the symmetry of Hq, the 
two-sideband state is symmetric [T] and can be written 
as qq = D 2 {a)S 2 { 0 iy^n{N)^iy_Q{N)Sl{^)Dl{a), where 
D 2 {a) = exp{[cr(a^j^ + a^j^)—h.c.]/v/2} is the symmetric 

displacement operator and 5'2(<a) = exp(^a^j^a^^ — h.c.) 
the two mode squeezing operator and (N) is the ther¬ 
mal state of mode a±Q with N average photons [2]. The 
state qq belongs to the so-called class of the two-mode 
squeezed thermal states, generated by the application of 
5'2(0^2(^) thermal states with (in general) dif¬ 

ferent energies. In order to test our experimental setup, 
we acted on the OPO pump and on the phase modulation 
to generate and characterize three classes of states: the 
coherent (a % 0 and = 0), the squeezed (^, % 0 

and Of = 0) and the squeezed-coherent (a, % 0) two¬ 

mode sideband state. We now consider the mode opera¬ 
tors: 


6(t, 0|f2) = Us, and 6(t, 7r/2|f]) = Ua, (2) 

which correspond to the symmetric {S) and antisym¬ 
metric (A) combination of the sideband modes, re¬ 
spectively, and the corresponding quadrature operators 
Qk = pk = and = 

X±7r/4(L ^/c|^): = a, s, with Ts = 0 and Ta = 7r/2. 

In the S/A modal basis, the first moment vector of qq^ 
reads R! = ((^s), (Ps), (^a), {Pa})^ and its 4 x 4 GM can 
be written in the following block-matrix form: 



where [38] : 

is the GM of the mode /c = a, s, e/ = (/g/a) “ (4)(G)5 
5ii = (/gTa) ~ (4)(L) with /,r = q^p and I % 1. The 
matrix elements of <Jk can be directly measured from the 
homodyne traces of corresponding mode |39|, whereas 
the entries of cr^ cannot. However, the information about 
ei can be retrieved by changing the value of the mixer 
phase to T = ±7r/4. In fact, it easy to show that [38l 
EOl m] e; = 5 ((^i) - (14) - I = Q^P, where 

q± = Xo{t, ±7r/4|f2) and p± = X„/ 2 {t, ±7r/4|r2). 

We now focus on dij. Given the state qq, but with 
different thermal contributions, these elements are equal 
to the energy unbalance between the sidebands (with¬ 
out the contribution due to the displacement that does 
not affect the GM) [39|, namely, 6qp = —6pq = ANq = 
(A^+o — A^-o), which cannot be directly accessed by the 
spectral homodyne detection alone. To overcome this is¬ 
sue, a resonator detection method has been proposed and 
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/ 1.00 ± 0.03 0.1 ± 0.1 0.04 ± 0.09 ( 0.1 ± 0 . 1 ) x 10 -^ \ 

1.03 ± 0.3 -( 0.1 ± 0 . 1 ) X 10-4 0.0 ± 0.2 

“ 1.01 ± 0.01 0.0 ± 0.02 
V 1.00 ± 0.02 

FIG. 2: Homodyne traces referring to the coherent two¬ 

mode sideband state and the reconstructed R' and cr'. The 
purities of the modes S and A are /ig = 0.99 ± 0.02 and 
/ia = 0.99 di 0.01, respectively. Only the relevant elements are 
shown. 


demonstrated in Refs. [Dill]. In our case we can exploit 
the error signal from the PDH stabilization to check the 
symmetry of the sideband state and also to measure the 
presence of some energy unbalance of the two sidebands, 
leading to non-vanishing 6if. More in details, given the 
cavity bandwidth, the PDH error signal allows to mea¬ 
sure the unbalance as [39] ANq = — where 

r±Q are the relative transmission coefficients associated 
with the two sideband modes and Nq = can 

be obtained from the (reconstructed) diagonal elements 
of (Tg and (Ta [39l|42]. 

Experimental results - Given the state po, the full 
reconstruction of the CM requires the measurement of 
the quadratures of modes as,aa, and a± = 6(t, ±7r/4|0). 
Once the mode has been selected by choosing the suit¬ 
able mixer phase 4^, the LO phase 0 was scanned from 0 
to 27r to acquire the corresponding homodyne trace.The 
statistical analysis of each trace allows to reconstruct the 
expectation value of the moments of the quadrature re¬ 
quired to reconstruct the CM a' and the first moments 
vector R'. 

Figures and [^ show the experimental spectral ho¬ 
modyne traces corresponding to the coherent, squeezed 
and squeezed-coherent two-mode sideband states, respec¬ 
tively. In the same figures we report the correspond¬ 
ing cr' and R'. All the reconstructed cr' satisfy the 
physical condition cr' ift > 0 where ft = icTy 0 cr^, 
(Ty being the Pauli matrix [2]. This implies that the 
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\ - 0.01 ± 0.06 / 


/ 0.50 ± 0.03 0.03 ± 0.06 0.02 ± 0.08 ( 0.1 ± 0 . 1 ) x 10 -^ \ 

4.2 ± 0.2 -( 0.1 ± 0 . 1 ) X 10-2 0.02 ± 0.07 

^ “ 0.49 ± 0.03 0.03 ± 0.05 

V 4.2 ± 0.2 

FIG. 3: Homodyne traces referring to the squeezed two¬ 
mode sideband state and the reconstructed R' and cr'. The 
noise reduction is 3.1 ± 0.3 dB for both the modes S and A, 
wheres their purities are /Xg = 0.6800.07 and /Xa = 0.6700.02, 
respectively. Only the relevant elements are shown. 


• Two-mode eoherent state 



R=(0.05±0.06, 2.18 + 0.05, 0.01 ±0.06, 2.23 + 0.057 


M .00+ 0.02 

0.0 + 0.1 

0.00 + 0.02 

0.1+0.1 ^ 



0.0 + 0.1 

1.02 + 0.02 

0.0 + 0.1 

0.01 +0.02 



0.00 + 0.02 

0.0 + 0.1 

1.00 + 0.02 

0.0 + 0.1 



, 0.1+0.1 

0.01 +0.02 

0.0 + 0.1 

1.02 + 0.02, 



• Two-mode squeezed state: 


R = (0.02 ± 0.04, 0.03 ± 0.04, 0.03 ± 0.04, 0.05 ± 0.04^ 



f 2.3 + 0.1 

0.00 + 0.06 

-1.8 + 0.1 

0.05 + 0.06^ 



0.00 + 0.06 

2.3 + 0.1 

0.01 +0.06 

1.8 + 0.1 



-1.8 + 0.1 

0.01 +0.06 

2.3 + 0.1 

0.00 + 0.06 



.0.05 + 0.06 

1.8 + 0.1 

0.00 + 0.06 

2.3 + 0.1 , 



• Two-mode squeezed-eoherent state: 


R = (-0.09 ± 0.05, 4.02 ± 0.06, -0.03 ± 0.05, 4.02 ± 0.06/ 



1 2.4 ±0.1 

0.1 ±0.2 

-1.8±0.1 

0.0 ±0.2 ^ 



0.1 ±0.2 

2.3 ±0.1 

-0.2 ±0.2 

1.8±0.1 


o-n = 

-1.8±0.1 

-0.2 ±0.2 

2.4 ±0.1 

-0.1 ±0.2 



. 0.0 ±0.2 

1.8±0.1 

-0.1 ±0.2 

2.3 ±0.1 , 



TABLE I: Reconstructed first moment vectors R and CMs 
ctq of the two-mode sideband states gn corresponding to the 
states of Figs. EE and[^ respectively. 


modes S and A represent the same local quantum state, 
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\ 0.04 ±0.06 


/ 0.53 ± 0.04 -0.1 ±0.2 0.02±0.03 (0.1 ± 0.1) x 10-^ \ 

4.1 ±0.2 -(0.1 ±0.1) X 10-2 0.2 ±0.4 

^ “ 0.57 ±0.03 -0.09 ±0.05 

V 4.2 ±0.2 J 

FIG. 4: Homodyne traces referring to the squeezed-coherent 
two-mode sideband state and the reconstructed R' and a'. 
The noise reduction is 2.7 ± 0.3 dB for the S mode and 2.4 ± 
0.2 dB for the S mode, wheres the purities are /is = 0.68±0.07 
and /ia = 0.64±0.02, respectively. Only the relevant elements 
are shown. 


namely, erg = CTai this is in agreement with our mea¬ 
surement within statistical errors, as one can check from 
the Figs. Hi and|^ Furthermore, the diagonal elements 
of the off-diagonal blocks are zero within their statistical 
errors, in agreement with the expectation for a factorized 
state of the two modes. 

We should now calculate the corresponding CMs in 
the modal basis d+o and d_o of the upper and lower 
sideband, respectively. Because of Eqs. we can write 
aQ = a' S and R = R\ where 


is the symplectic transformation associated with the 



s = ^( “ “ 

V2 \ —'iCTy icr. 


mode transformations of Eqs. <§• The results are sum¬ 
marized in Table m Whereas the reconstructed two-mode 
sideband coherent state is indeed a product of two coher¬ 
ent states, the other two reconstructed states exhibit non- 
classical features. In particular, the minimum symplec¬ 
tic eigenvalues of the corresponding partially transposed 
CMs m [44] read A = 0.50 ± 0.02 and A = 0.55 ± 0.03 
for the two-mode squeezed and squeezed-coherent state, 
respectively: since in both the cases A < 1, we conclude 
that the sideband modes are entangled. 

Concluding remarks - In conclusion, we have presented 
a measurement scheme to fully reconstruct the class of 
symmetric two-mode squeezed thermal states of spectral 
sideband modes, a class of states with Gaussian Wigner 
functions widely exploited in continuous variable quan¬ 
tum technology. The scheme is based on a homodyne de¬ 
tection and active stabilization, which guarantees phase 
coherence in every step of the experiment, and on a suit¬ 
able analysis of the detected photocurrents. We have 
shown that by properly choosing the electronic mixer 
phase it is possible to select four different combinations 
of the upper and lower sideband which, together with the 
information form the PDH error signal, allows to recon¬ 
struct the elements of the covariance matrix of the state 
under consideration. The scheme has been successfully 
demonstrated to reconstruct both factorized and entan¬ 
gled sideband states. 

In our implementation we have used two electronic 
mixers and retrieved information about two modes at 
a time. It is also possible to use four mixers and extract 
information about the four modes at the same time. The 
method is based on a single homodyne detector and does 
not involve elements outside the main detection tools of 
continuous variable optical systems. As such, our pro¬ 
cedure is indeed a versatile diagnostic tool, suitable to 
be embedded in quantum information experiments with 
continuous variable systems in the spectral domain. 
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SUPPLEMENTAL MATERIAL 

1. Covariance matrix elements of the two-mode squeezed thermal state 


In this section we explicitly show how we can calculate the elements of the CM cr', given in Eq. (3) of the main 
text. We recall that, according to our definitions: 

^ln^ a+Q(t) + a-n(t) . a+ct{t) - a-ct{t) . 

b[t,0\il) = - — - = ds, and b[t,7r/2\il) = t - — -= (6) 

V 2 V 2 

therefore the quadrature operator Xo{t^ 4^|Q) = 5(t, ^|^^) -T 5’*’(t, ^|^^) e*^ can be written as: 

Xo(t, 4^|Q) = cos A/ [qs cosO Ps sin6>] + sin 4^ [g^ cosO Pa sin6>] . (7) 

If we set 4^ = 0, we have: 


Xoit, 0|e) = (Zs = as + at = ^ {q^,} - {qsf, 


X^/2{t,0\Q} =Ps = i{al - as) = ^ ^ (Ps) “ (Ps)t 

v±^/4(i,o|e) = => ^{qsPs + PsQs) - {qs){ps), 


(8b) 
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for ^ = 7r/2 we obtain: 


Xo(i,7r/2|ll) = q^ = a^ + al = ^ (ql) - {q^f, (9a) 

W/2(i, 7r/2|fJ) = Pa = *(4 - fla) = ^ {pD “ (Pa)t (9b) 

^±7r/4(i,7r/2|fi) = ^ ^(^aPa +Paga) “ (9a)(Pa), (9c) 

On the other hand, if we set ^ = ±7r/4 we find: 

Xo(t,±7r/4|n) = ^5^, and X,/ 2 (i, ±7r/4|n) = (10) 

ad we have the following identities: 

{Xl{t,n/m - Xl{t, -n/m) = 2 ( 9 a 9 s) = e„ ( 11 ) 

7r/4|n) - Xl/^{t, -7r/4|ll)) = 2(paPs) = tp- (12) 


As mentioned in the main text, it is not possible to calculate the elements and 5pq directly from the spec¬ 
tral homodyne traces [T]. However, when the state under consideration is a two-mode squeezed thermal state 
Qn = D 2 {a)S 2 {C)v+a{Ni) ® v-Q,{N 2 )Sl{C)D 2 {a), where D 2 {a) = exp{[Q;(a(^f 2 + aLn) - h.c .]/\/2} is the symmet¬ 
ric displacement operator and S 2 {oi) = exp(^a^Qa^Q — h.c.) the two mode squeezing operator and is the 

thermal state of mode a±Q with N average photons [2], we can calculate 5qp and 5pq as follows. 

Since the covariance matrix does not depend on the displacement operator, we can assume a. Furthermore, it is 
useful to introduce the following parameterization: we define the squeezed photons per mode A^sq = sinh^ r, the total 
number of thermal photons A^th = A/^i + A ^2 5 and the thermal-photon fraction i?th = A'l/A^th- Thereafter, the energies 
of the two sidebands are given by: 


A/'+O — A^sq(l + Ath) + ^thAth, 

A_o = Asq(l -h Ath) + (1 — ^th)Ath, 


respectively, and thus: 


A+o + A_q — 2A^sq + Ath(l + 2A^sq) and A+f^ — A_f^ — A^th(2i?th “ !)• 


The covariance matrix associated with have the following block-matrix form: 

( CcT,\ 

(Tn= [ 

\CcT, BI J 

where az is the Pauli matrix and: 

A = 1 + 2A^sq(l + Ath) + 2i?thAsq, 

H = 1 + 2A^sq(l + Ath) + 2(1 — i?th)Asq, 

C = 2(1 + iVth)t/A^sq(l + iVsq) . 

The corresponding measured covariance matrix reads (see the main text for details): 

_ n{A + B)I + Ca, {N+n-N.n)i(Ty \ 

\ {N+a-N_a)i(Ty liA + B)I + CaJ ’ 


(13) 

(14) 


(15) 


(16) 


(17a) 

(17b) 

(17c) 


( 18 ) 


ay being the Pauli matrix. Note that while a' is always symmetric, a^ can be also asymmetric. 
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2. Retrieving the energy unbalance of the two-mode squeezed thermal state 

In order to determine the energy difference between the sidebands, we should measure the cavity bandwidth and, 
by exploiting the error signal of the PDH [3], we can assess the relative cavity transmission coefficients t±q = 
7±o/(7+o + T-n) associated with the two sideband modes, where T±q are the actual transmission coefficients. 
Therefore, the energy difference can be obtained as: 


= + (19) 

In general, given the covariance matrix tr of a Gaussian state, the total energy can be obtained from the sum of its 
diagonal elements [cr]kk as (without loss of generality we are still assuming the absence of the displacement): 

kk - 'i- ( 20 ) 

^ k=l 

Experimentally, we can find the total energy from the first and second moments of the operators in 

Eqs. ^ and in Eqs. ([^, which are measured from the homodyne detection. 
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